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Abstract — Symmetry is a fundamentally important concept 
in many branches of physics. In this work, we discuss two 
types of symmetries, external symmetry and internal symmetry, 
which appear frequently in controlled quantum spin chains and 
apply them to study various controllability problems. For spin 
chains under single local end control when external symmetries 
exists, we can rigorously prove that the system is controllable 
in each of the invariant subspaces for both XXZ and XYZ 
chains, but not for XX or Ising chains. Such results have direct 
applications in controlling antiferromagnetic Heisenberg chains 
when the dynamics is naturally confined in the largest excitation 
subspace. We also address the theoretically important question of 
minimal control resources to achieve full controllability over the 
entire spin chain space. In the process we establish a systematic 
way of evaluating the dynamical Lie algebras and using known 
symmetries to help identify the dynamical Lie algebra. 

Index Terms — quantum control, spin chains, symmetry, sub- 
space controllability 

I. Introduction 

Controllability is a fundamental concept in control theory 
in general, and control of quantum systems in particular. Any 
quantum system with a sufficient number of controls becomes 
fully controllable l[l)-l|4). Therefore we are most interested in 
the problems where the system has only a limited number 
of controls and often limited controllability (e.g. subspace 
controllability). Such limited controllability is usually due to 
the existence of symmetries in the Hamiltonians Q, |j6), which 
restrict the dynamical Lie algebra (DLA) of the system fT^. 
Previous literature on quantum controllability has mainly 
focussed on the cases where either the system is fully control- 
lable (hence implying universal quantum computation |8|), or 
not fully controllable but with a DLA that scales linearly or 
quadratically with the system size. In contrast, in this work, we 
would like to study systems that are not fully controllable but 
with a DLA large enough for universal quantum computation. 

There are simple criteria for controllability of bilinear 
systems in terms of the Lie algebra rank condition |9| similar 
to the Kalman rank condition for linear systems. However, 
verifying controllability for quantum systems is challenging. 



not least because the dimension of the DLA associated with 
a multi-partite quantum system usually grows exponentially 
in the number of particles (such as qubits). This exponential 
scaling makes it impossible in most cases to verify the Lie 
algebra rank condition numerically. It is therefore important to 
have general algebraic controllability results for certain classes 
of systems such as spin chains with a few controls of a certain 
type. In this paper we derive such results for spin chains 
with isotropic and even more anisotropic couplings. Unlike 
spin chains with Ising-type coupling such systems are usually 
controllable with very few controls acting on a small subset 
of spins. However, controllability is limited by the existence 
of symmetries in the Hamiltonians. For instance, it has been 
shown using the propagation property that Heisenberg chains 
are fully controllable given two non-commuting control acting 
on the first spin |1| but not when there is only single control 
acting on the first spin. In the latter case the controlled system 
has symmetries and decomposes into invariant subspaces |6|, 
preventing full controllability. However, it has been observed 
that such systems appear to be controllable on each invariant 
subspace, in particular, the largest excitation subspace, whose 
dimension scales exponentially with system size |10|. 

In this work we give a rigorous proof of this subspace 
controllability result for XXZ chains and then apply similar 
techniques to discuss the subspace controllability of a general 
XYZ spin chain. This system is interesting as it provides 
arguably the simplest model of a universal quantum computer 
one could imagine: a physical Hamiltonian with a single 
control switch to do the computation. We further show that 
the same result does not hold for XX chains, where a single 
control acting on the end spin in a chain can only give 
controllability on a subspace whose dimension does not scale 
exponentially with system size. In this case additional controls 
are needed, and we discuss the minimal local control resources 
for full controllability in this context. 

This paper is organized as follows: in Section |Il] we intro- 
duce different types of spin chains and define two fundamental 
types of symmetry, external and internal, and their relations to 
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controllability. In Section|lIl] we present a complete discussion 
on spin chains under a single end control, and rigorously 
prove that for both XXZ and XYZ chains, the system is 
controllable in each invariant subspace, and that this result 
is robust if the control field has a leakage on the neighboring 



spins. In Section IV we investigate the XXZ or XYZ chains 
for various types of two controls and we find the minimal 
control resources for full controllability on the entire Hilbert 
space. In Section [V] we study the dynamical Lie algebra for 
an XX chain subject to a single end control and investigate 
the controllability for an XX chain subject to two and three 
controls. 

II. Model and Basics 

For a quantum system composed of N spins, we denote 
the standard Pauli operators by X, Y, Z and the local operator 
R acting on the k-th spin by Rk, i.e., R^ = I ■ ■ ■ IRI ■ ■ ■ I, 
where / is the identity on a single spin. 

System Hamiltonian: We consider a spin network com- 
posed of N spin-i particles with spin-spin interaction charac- 
terized by the following Hamiltonian 
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with the special cases am„ = 6^ 
isotropic Heisenberg coupling, a,. 



^mn Zra Ztl ( 1 ) 

I corresponding to 
to XXZ-coupling, 



to XX-coupling, and a.,„„ — = 
to Ising coupling. For XXZ-networks it is convenient to set 
7mn = a-mn = and c„„ = '«7mn- We also require all 
couplings (to, n) in ([T]l form a connected graph. 

The constants amm^mmCran determine the coupling 
strengths between nodes to and n in the network. Special 
cases of interest are chains with nearest-neighbor coupling, 
corresponding, e.g., to linear qubit registers in quantum in- 
formation processing, for which a„i„ = 6„i„ = c„i„ = 
except when to = rt ± 1. A network is uniform if all non- 
zero couplings are equal, i.e., e {0, a}, 6,„„ e {0, 5} and 
Cmn € {0,c}. Every spin network has an associated simple 
graph representation with vertices {!,.... iV} determined by 
the spins and edges by non-zeros couplings, i.e., there is an 
edge connecting nodes to and n exactly if 7m„ ^ 0. 

Controllability: The controlled quantum dynamics we are 
interested in is characterized by the following Schrodinger 
equation: 



P = 



(2) 



where Hq is the system Hamiltonian in ^ and Hj, j — 
1,2, ••• ,TO is a series of control Hamiltonians with time- 
varying amplitudes fj{t). We define the system to be con- 
trollable if the dynamical Lie algebra C generated by iHj, 
j = 0, 1,-- - ,TO is equal to the largest Lie algebra u(2^) 
or 5u(2^). The definition of controllability is very intuitive: 
it can be shown that if the system is controllable, then any 
unitary process U G SU(2^) can be generated from ^ under 
certain control sequence in finite time; if £ C su(2 ), then 
there exists some unitary gate U G SU(2^) that can never 



be generated under Q ||7|. The concepts of controllability 
and dynamical Lie algebra are very important for both theory 
and control applications, as they characterize the reachable 
set of the control dynamics and have answered the question 
whether a given control task can be achieved or not. However, 
calculating the dynamical Lie algebra can become extremely 
difficult or even intractable as N increases. Therefore, we hope 
to use other properties of the Hamiltonians to infer information 
about controllability, and symmetry does play such a role. 

Symmetries: We consider two types of Hamiltonian sym- 
metries: external symmetry and internal symmetry [3|. 

Definition 1. Let Hj, j = 0,1, ■ ■ ■ , m, be a set of Hamilto- 
nians for a given quantum system. If there exists a Hermitian 
operator S such that [Hj , S*] = for all j then S is called 
an external symmetry for the Hamiltonians; assuming Hj 
are trace-zero, if there exists a symmetric or antisymmetric 
operator S such that Hj S + SHj ~ for all j, where HJ is 
the transpose of Hj then S is called an internal symmetry. 

From the definition, external symmetry implies that all Hj 
can be simultaneously diagonalized, while internal symmetry 
implies that the dynamical Lie algebra £ generated by Hj is 
a subalgebra of the orthogonal algebra so (2^) or symplectic 
algebra sp(2^) |1 1) . In both symmetry cases, £ is strictly 
smaller than 5u(?'^) and the system is not controllable. It 
is useful to investigate which operators can be the external 
symmetries. 

Example 1. For the system Hamiltonian ([T]), a simple class 
of symmetry operators is of the form S = AiA2...An, 
where A^. is a local operator on the fc-th spin, i.e., A^ ~ 
aiX + a2Y + a^Z + a^I . [iJoi'S'] — then requires 
[XrnXn, ArnAn] = for any connected link {m,n) in ([T]l, 
which shows that the nontrivial external symmetry operators 
are X1X2 ■ ■ ■ Xn, Y1Y2 ■ ■ ■ Iat and Z1Z2 ■ ■ ■ Zj^, which are 
often known as the parity symmetry. Hence, if the control 
Hamiltonians only contain local Pauli operators in one direc- 
tion, such as Z direction, with Hj = Fj{Zi, . . . , Zn), then 
Sp = Z1Z2 ■ ■ ■ Zm is the corresponding parity symmetry and 
all Hamiltonians are invariant in each of the two eigenspaces 
of Sp with parity +1 and —1. 

Example 2. If the system Hamiltonian ([T]i is of XXZ type 
then defining Se = (^j + -^)/2, we have [Ho,Se] = 0. 
Physically, Se represent the total number of excitations, and 
has + 1 distinct eigenvalues, ranging from n = to 
n = N, corresponding to different numbers of excitations 
in the network. If the control Hamiltonians only contain Z 
operators, i.e., Hj = Fj{Zi, . . . , Z^) then 5*6 defines an 
external symmetry, called the excitation symmetry, and all 
Hamiltonians are block-diagonalized on the iV + 1 invariant 
subspaces, as illustrated in Fig. [T] 

Example 3. A non-identity element 11 in the permutation 
group defines a permutation symmetry of the spin network 
if all Hamiltonians are invariant under the permutation 11 of 
the spin indices [6|. In particular, for a single-local-control 
problem such as Hi = Zk on the fc-th spin, permutation 
symmetry means that the index k much be fixed under the 
permutation 11 (Fig.|2]i. In fact, 11 induces a symmetry operator 
5n which commutes with both the system and the control 
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control Hi = Z\ on the end spin of a XXZ chain, with the 
symmetry operator S^. as the total excitations. 



d = 3 



Fig. 1: For = 4, the XXZ network Hamiltonian Rq and 
local control iJi = are simultaneously block-diagonalized 
in 5 excitation subspaces, "Kd, d = 0, . . . , 4. 



^3 




Fig. 2: Different configurations ||5J of controlled spin networks 
with a local control li\ = Zk on the k-\h spin indicated 
by the arrows (the indexed circles represent different spins, 
with edges as their couplings): (a),(b),(c) have permutation 
symmetry, (d) does not. 



Hamiltonians and hence defines an external symmetry. More- 
over, since Sn commute with Se as defined in Example 2, 
also induces external symmetries on each excitation subspace 
of Se, i.e., all Hamiltonians can be further block-diagonalized 
in the excitation subspaces. 

Having found all external symmetry operators of the Hamil- 
tonians, the entire Hilbert space can be decomposed into 
T-L — ®d=i^d, where quantum dynamics is invariant on 
each T-Ld, which cannot be further decomposed. The asso- 
ciated dynamical Lie algebra L must be a subalgebra of 
©£^j^u(dimH(j). Although the system is not controllable on 
the entire space, it may still be controllable on each T-Ld- In 
the following, we show that this is indeed true for single local 



III. Single Local End Control 

One of the simple but important configurations of a spin 
network is a spin chain, which is the main subject of the paper. 
We first consider a spin chain with a single local control at 
the end of the chain. Without loss of generality, we assume 
the control field is in Z direction. The corresponding control- 
lability result depends on whether the spin-spin interaction on 
the other two directions are equal or not, i.e., whether the spin 
chain is of (I) XXZ type or (2) anisotropic XYZ type. 

A. XXZ Chain 

For an XXZ chain with spin number TV, under the end 
control in Z direction, the system and the control Hamiltonians 
are written as: 

N 

Ho = J2 (^^■^^■+1 + ^^■^^■+1 + ^jZ,Z,+i) (3a) 
i 

Hi = Zi (3b) 

As discussed in previous section, the excitation operator 
Se = {Zj + /)/2 is an external symmetry, and the entire 
Hilbert space is decomposed into % = ©fcLo^fe Hk as 
the invariant subspace with k excitations, i.e., it is generated by 
the computational basis vectors with k number of I's, where 
the two single-spin basis vectors are denoted as |0) and 
Hence dim(Hfe) = (^) = dAr.fc. For example, for = 4, 
H2 is expanded by |0011), |0101), |0110), |1001), |1010) and 
1 1100) with dim('H2) = ^4,2 = 6. Due to Se, the controlled 
system (|3]l is not fully controllable on the whole space, but 
it is controllable on each Hk- As an application, for \j < 
when Hq represents an antiferromagnetic chain, and we can 
easily prepare the system into the ground state ipo, which is 
in the largest excitation subspace 'H[n/2\ at t = to- Then, 
by applying a single control Zi with amplitude f{t) derived 
from optimization, we can generate the total Hamiltonian 
H = Ho + f{t)Zi to drive the system into an arbitrary target 
state in T-LyN/2\ at a later time t — tp-^Vi particular, we can 
generate perfect entangled pairs between the two end spins of 
the chain, which is an important quantum resource for many 
applications such as quantum communication or measurement- 
based quantum computing fTOl. 

Next we rigorously prove that under the control dynamics 
with the Hamiltonians in Q, the system is controllable in each 
T-Lk, and particularly in T-LyN/2\ - By definition of controllabil- 
ity, it is sufficient to show that Hq and Hi generate u{d]si^k) 
on each T-L^. Since T-L = ©^g'Hfe, the associated dynamical 

Lie algebra C C Ct = ®j^o^('^^-fc)- ^'^^^ °f '■^^ proof 
is to determine all independent operators generated in C and 
then evaluate dim(£) in order to identify C. 

Since a Lie algebra is also a real vector space, we can drop 
some factors in the calculation and use linear combinations. 
We denote such (trivial) steps in the derivation by First of 
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all, we derive the following commutation relations: 

\X\X2 + Y\Yi, X1Y2 — Y1X2] Z2 — Z\ ^ Z2 



Continuing this process, we can generate all Zj, XjX^ + 
YjYk and XjYk — YjXk (with details in appendix |A]). For 
brevity purposes we will only focus on the XX + YY terms 
and not write down the XY — YX terms explicitly, since 
one operator can always be generated from the other An 
operator is called a fc-body operator if it contains k nontrivial 
factors, i.e., those comprised of X, Y 01 Z Pauli operators. 
For example, Z1Z3 is a 2-body operator, while (^1 — ^2)^3^4 
is a 3-body operator. Denoting Mk as the set of all fc-body 
operators in C, we list all elements in Mk and evaluate 
rank : 

(1) Afi = {Zk}; 

(2) M2 = {XjXk + YjYk, Z,Zk}; 

(3) M3 = {{XjXk + Y,Yk)Zrn, (Zj - Zk)Z^Zn}; 

(4) M4 - {{XjXk + Y,Yk)ZraZ^, 

{XjXk + YjYk){XmXn + YjnYn), {Zj — Zk)ZiZ„iZ„}. 

(5) M^ \^(^X^-^Xjji^ ~\- YjYi-^Yjyi^ Zjyi^Zyyi^Zrjy^r-^^ 
{X^-^X^l^ -\- Yyyi-^Yyyi^^(yX'fyi.^Xui^ Yyy^^Yyy^ Z yyi ^ ^ 
{,Zmi Zjji^^ Zjyi^ Zjyi^Zjyi^^ ZjyiQ^ , 



(£) Mf. when £ is even, we can generate: 



i^X'fji-^X'fYi^ Yrfyi-^Yrfyi^^ • • • (^Xjji^ — iXjji^ -\- Y^^^iY^^^ 

{Zjji-^ Z yyi^^ Z jyi.-^ ' ' ' Z 

when £ is odd, we can generate: 

(.Xjji-^^Xjyi^ -\- YfY^-^Yyy^^^Zj^-i^ ' ' ' Z ^ 

(.Xjyi-^Xjji^ -j- Yjyi-^Yjyi^^(^Xjyi.-^X^^ "t" Yyy^^Y^fi Z yy^ ^ ' ' ' Z yfi ^ 



However, not all of the elements in i?jv,p are linearly 
independent. For example, for = 4 and p = 2, we find 

[x^X2 + rir2)(X3A4 + Y^Y^) - (X1X3 + Yiys) 

(X2X4 + ^2^4) - (Xin - YrX^){X2Y^ - 1^2X3), 
{X^Y2 - Y^X2){X^Y^ - ^3X4) - (X^Y^ - FiXg) 
(^2^4 - y^X^) = (X,Y4 - YiX^){X2Y3 - Y2X3), 
{X1X2 + YiY2){X3Y4 - Y3X4) - (X1X3 + YiYs) 
{X2Y4 - Y2Xi) = {X1X4 + YiY^){X2Y3 - Y2X3) 

Similarly we can write down the other dependence relations. 
Altogether there are only 1/2! of all 2-pair operators that are 
linearly independent. In general, we will prove that only l/p\ 
of all p-pair operators are linearly independent, and 

'A^\ /N - p 

,-P/ V P 

However, directly proving Q is very difficult as the linear de- 
pendence relations can become very complicated for large A^ 
and p. Fortunately, we can convert this problem to evaluating 
the rank of a set of polynomials on complex field C (with 
details in appendix IbI). Therefore, for £ > 2, 



i^XjYi-^X'^^ -j- Yfyi-^Yfyi^^ • • • (A^^^^_2 A^m^— 1 ~^ ^171^— 2^7711— 1) Zj^ 



ma ■ ■ ■ Zyfi^ 



Next, in order to get rank(Af£), we first evaluate the number 
of the operators in the form 

{Xm-^Xjji^ + Yjn-^Yjyi^] ■ ■ ■ (A'„j2p-i ^m2p ^" ■^n2p- 1 ■^^2p ) J 

which contains p pairs of {XX + YY) or {XY - YX) 
operators. We will call them p-pair operators. For a given 
N and p with A^ > 2p > 0, we denote the set of p- 
pair operators as Ejy p. For example, for A^ 2p = 6, 
{X1X2 + YiY2){X3Y^ - Y3Xi){X5Xe + Y^Ye) is a 3-pair 
operator in Eq ^. Then the size of the set E^^p is obtained by 
simple combinatorics as 

2PfN\fN-2\ fN-2{p-l)\ fN\fN-p 



rank(i?jv,p) 



(4) 



rank(Af£) = rank(i?Ar^p) 
p=i 



£-2p 



1. 



In the above, the extra combinatorial factor arises from differ- 
ent choices for putting the Z terms. After some simplification, 
for both AT = 2to and A^ = 2m + 1, 



dim£ =^rank(M£) 

e=i 

=:Tank{EN,i)2^-^ 
+ 2^ - A^ + 1 

lN/21 



- rank(£'Ar,„)2 



N-2m 



j^\2N-2p 



pP{N-2p)l 



E 



f2N 



N + 1 



2N 

N 



where the last equation is shown in Lemma |2] in the appendix. 
As discussed earlier, C C Ct = ®k=o'^i'^N,k)^ with 

N N 

dim(/:T) = X! dim(u((iAr,fe)) = ^ 

4=0 2=0 

All fc-body Z-type operators, fc = 1, • • • , A^, generate a Cartan 
subalgebra C in Ct, with dim(C) = 2^. Notice that since 
we can only generate coupled fc-body Z-type operators, such 
as {Zrm - Zjn2)Zm3 ■ ■ ■ Z,nt 1" C, the rank of all Z-type 
operators in £ is 2^ — A^ + 1, i.e., there are A^ — 1 independent 
Z-type operators not included in C, but included in Ct- We 
hence have: 

/2Ar^ 



dim(£) < dim{CT) - N + 1 



N + 1 = dim(£) 



It means that dim(£) achieves the allowed maximal value, 
which is true only when C is isomorphic to u{d]y k) on each 
T-Lk- Hence, we have proved the following theorem: 

Theorem 1. For an XKZ chain of length N with a single 
local control on the end spin in Z direction, the system 
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is controllable on each 
subspaces. 



of the N + 1 invariant excitation 



In particular, this theorem holds for anti-ferromagnetic 
Heisenberg chains, which rigorously justifies the numerical 
findings in |10|. Moreover, as 'HyN/2\ is exponentially large 
as N increases, it can be used as a resource for universal 
quantum computation. For instance, we can encode qubits as 
a|01) + /3|10), thereby performing universal quantum com- 
putation in T-LyN/2\- This is a remarkable observation: we 
have found a system where quantum computation can be 
achieved with a single switch, and where both the system 
and control Hamiltonian are physical, e.g. consist of nearest- 
neighbor two-body interactions, which are very common in 
physics. It provides possibly the simplest and most elegant way 
of achieving quantum computation so far (leaving efficiency 
issues beside |4|). Having only a single switch we avoid the 
experimental difficulty of quickly changing field directions. 



operators: 

(1) Mi: Zk-, 

(2) M2. PjPk, where P can be X or Y, and Z^Z^; 

(3) M3 



PjPkZn 



and Zu Zm Zn 



(£) Mf. when £ is even, we can generate: 

P P 7 7 

PmiP'rn2Pm^PraiZra^ ' ' ' ^ 



P P 



-iPrr 



When i is odd, we can generate 

■■z 



P P P P 



mi 



B. XYZ chain 

For XYZ chain under Zi control with 

N 

^0 = ^ flj^j^j+i + bjYjYj+i + CjZjZj+i 
j 

Hi=Z, 



(5a) 
(5b) 



where aj 7^ bj, does the subspace controllability still exist? 
As discussed in Example 1, there exists a parity symmetry 
Sp — Zi - ■ ■ Zff satisfying [Hq, Sp] — [Hi,Sp] — 0, with two 
invariant subspaces Hi and H-i, corresponding to eigenvalues 
±1 of Sp. We will show that the Hamiltonians cannot be 
further block-diagonalized on each of the two subspaces, 
and the system is controllable on each of them. Notice that, 
compared to the XXZ chain, the number of invariant subspaces 
for XYZ chain has reduced from TV + 1 to 2, which is 
not too surprising as we have broken the symmetry between 
X and Y directions from XXZ to XYZ type, and some 
symmetries should disappear. In the following we will identify 
all operators in £ generated by and Hi: 

[Zi,Ho] aiYiX2 - biXiY2 ^ 01X1X3 + 61^1^2 
[aiYiX2 - biXiY2,aiXiX2 + biYiY2] 
{al + bl)Zi-2aibiZ2 ^ Z2 

Continuing this process, we obtain all Zj, ajXjXj^i + 
bjYjYj+i' ajYjXj+i - bjXjYj+i and H2 = T,jCjZjZj+i. 
Next, we have 



and together with ajXjXj+i + bjYjYj+i we can decouple 
and get XjXj^i and YjYj^i. Similarly we can decouple 
and independently generate XjYj+i and YjXj+i. This a 
major difference from the XXZ case, where the XX and YY 
operators at neighboring locations cannot be decoupled. Due 
to such decoupling, we expect that the dynamical Lie algebra 
C generated by Hq and Hi will be larger than the XXZ case. 
Next, repeating the same generation process by calculating the 
commutators, we get the following set series Mk of fc-body 



Pmi Pm2 ' ' ' P'm^~2Pmi — lZrn^ 
Zfni Z^Yii ' ' ' Zrni 

Compared with XXZ chain, where we can only generate the 
coupled Z-type operator, such as {Zi — ^2)^3^4, for XYZ 
chain, we can separately generate Z1Z3Z4 and Z2Z3Z4. 
is be divided into two subsets: the set of P — Z-type operators 
and the set of Z-type operators, where each P — Z operator 
can contain 2p number of P's and N — 2p number of Z's, 
p — 1,. . . , [i/2\. Hence, following some basic combinatorics 
argument, we have: 



rank(M£) = ^ 



2''^(^^^]m\rp 



N -2p 
£~2p 



and the dimension of C{N > 2): 

L7V/2J 



k=0 



-2 = 22^^-1 - 2 



where we have used the identity 



LJV/2J 

E 

k=0 



LW/2J 

E 

k=0 



N 
2k + I 



= 2 



Since Hq and Hi are simultaneously block-diagonalized on 
Hi ® H-i, C must be a subalgebra of Ct = u(2^^i) ® 
u(2^^^). Moreover, since the fc-body Z operators in C are 
generated from the (fc + l)-body P — Z operators, C does not 
include two Z-type operators, the identity / and Sp, which 
are however included in Ct- Hence, we have 



dim(£) < dim(£T) -2 = 2 



2JV-1 



2 = dim(£) 



Hence dim(£) achieves the allowed maximal value, which 
is only true when restricted on each of the subspace Hi and 
H^i, C = u(2^~i) or su(2^-i). Noticing that Hq and Hi 
are trace-zero on Hi and H-i for > 2, we must have 
C ^ su(2^~i) on both Hi and H-i for > 2. When A^ = 2, 
it is easy to check that C — u(2) on Hi and H-i- Thus, we 
have proved the following theorem: 
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Theorem 2. For an XYZ chain of length N with a single 
local control on the end spin in Z direction, the system is 
controllable on each of the two invariant subspaces Hi and 

n-i. 

C. When Control Has a Leakage on Neighboring Spins 

The previous assumption of control on a single spin only 
holds in theory. In practice, it is difficult to apply a control field 
that only acts on a single spin without affecting its neighbors. 
Hence, a more realistic assumption is that the local end control 
has a leakage on the neighboring spins with Hi = 7i^i- 
We consider two common types of leakage: linear jj — —aj + 
j3 and exponential 7^ = e^^'^^^-* decays. In the following, 
we show that the subspace controllability results discussed so 
far are robust against such control leakage, i.e., when single 
control field Zi has some leakage on the neighboring spins, the 
system is still controllable in the invariant subspaces. Under 
the leakage assumption. 



N 



Hq - ^ ttjXjXj+i + bjYjYj+i + CjZjZj+i 



(6a) 



to the neighboring spins, the system is controllable on each 
of the invariant subspaces. 



IV. Minimal Controls for Full Controllability 

In previous section, we have provided a complete discussion 
of the control problem of spin chains with the least control 
degree of freedom, i.e., a single local control at the end of the 
spin chain. In general, as the number of controls increases, 
existing symmetries will disappear and the system will become 
fully controllable on the entire Hilbert space under a sufficient 
number of independent controls. Therefore, another interesting 
theoretical question is to ask when such transition happens 
from an uncontrollable system to a fully controllable one. 
Alternatively, we can ask what are the minimal controls that 
can make the chain fully controllable, which is the main 
topic of this section. We will base on the results in previous 
discussions and add more controls to the control systems under 
^ or 



Hi = ^7jZj, 

Defining adjoint action of Hi on Hi as Ad^i (^^o) = [^1 
and Aj — ajXjXj^i + bjYjYj+i. 



(6b) 



A. Controlling Zi and Xi 

In ||T|, it was proved by the propagation property that an 
XXZ chain with two independent controls Hi = Zi and H2 = 
ry^A]^i is fully controllable on the entire space. We can rederive 
^ this result from our analysis in previous section: observing 
'the operators generated by Hq and Hi, and writing down the 
operators generated by Hq and H2, it is easy to see that we can 

Ad^;^(i/o) - (71 - 72)'" Ai + • • • + (7fe_i - ^k^'Ak-i + jI'S^T'^^^ ^" '^"^"'^y ^^""^^ operators, k = l,...,Nm u(2^). 
^ Hence the system is fully controllable. 



Adgj(iJo) = (71 -72)^^1 + 
Ad(,^|(i/o)-(7i-72M + 

j(2fc)^ 
Hi 



we have 

+ (7/C-1 - 7/c)^ A-i 
+ (7/0-1 - lk)^A 



k-l 



where the coefficients in this expression can be denoted by 

f (71 - 72)^ ■ • • (7fe-i - 7fe)^ Ik \ 
(71 - 72)^^ ■ • • (7fc-i - 7fe)'* 



V = 



it 



V(7i - 72) 



2k 



(1) When the leakage of the local control is linear, i.e., 

7i ^ 7i+i = 7^ ^ 7^+1 for different j and i, we can generate 
the operator — a^X^Xk^i +6feyfcYfc+i from any two rows 
of V . Analogously, we can generate a^X^Ykj^i — b^Y^Xk^i 
and hence generate Zk~Zk+i- From Ak, we can also generate 
Bq = Ad^|(i/o) - Ik^k- From Zk - Zk+i and Bq, we can 
sequentially generate Aj and Zj, j — k — 1, k — 2, ■ ■ ■ ,1. 

(2) When the leakage of the local control decays nonlinearly. 



e-g-, 7i 



, we have - -fj+i ^ -fe - le+i ^ 7fc. 



and from the property of Vandermonde matrix, det(y) 7^ 0. 
Hence we can generate each Aj, j ^ 1, . . . , k. Together with 
Hi, we can decouple and generate Zj, j = 1, . . . , fc. 

Hence, in both cases, £ generated by Hq and Hi in (|6]l is 
the same as that generated by Hq and Hi = Zi. In general, 
for other types of nonlinear leakage, the above reasoning is 
valid for almost all choices of jj. Thus we have: 

Theorem 3. For an XXZ or XYZ chain of length N, under a 
single local control on the end spin in Z-direction with leakage 



Theorem 4. For an XXZ or XYZ chain of length N, with two 
local controls on the end spin. Hi — Zi and H2 — Xi, the 
system is controllable on the whole space. 



B. Controlling Zk and Xk 

In Theorem [4] we have shown that if we can fully control 
the end spin, then the system is controllable on the whole 
space. What if we can fully control one spin at other locations? 
We will prove that for a general XYZ chain two independent 
controls on the fcth spin in Z and X directions 



N 

Ho = ^ o-jXjXj+i + bjYjYj+i + CjZjZj+i 
i 

Hi = Zk, H2 = Xk 



are sufficient for controllability on the whole space, except 
when N ~ 2k + 1, where the Hamiltonians exhibit a mirror 
permutation symmetry with respect to the center kth spin. 
Specifically, for k < [-/V/2J, let us calculate the operators 
in C generated by the three Hamiltonians. 
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^{ak-iXk-iYk - bk-iYk-iXk) + {akYkXk+i - bkXkYk+i) 
-^{ak-iXk_iXk + bk-iYk-iYk) + {akXkXk+i + bi.YkYk+i) 
~^(«fe-i + bl-i)Zk - 2ak-ibk-iZk-i 
+ {al + bl)Zk - 2akbkZk+i 

+ 2{ak-iakXk-iXk+i + bk-ibkYk-iYk+i)Zk = Pk 
[Xki [Xk,Pk]\ — Pk ^ dk-iZk-i + rffe+i-^fe+i = Qi 

[Ql,Hz\ ■ ■ ■ ^ dk-2Zk-2 + dk+2Zk+2 = Q2 

Continuing this process, we can sequentially generate: 

dk-iZk-i + dk+iZk+i, dk-2Zk-2 + dk+2Zk+2, • ■ ■ , diZi + 
d2k-iZ2k-i, Z2k — Z2k-i, Z2k+i — Z2k, ■ ■ ■ ,Zn — Zn^i. as 
well as 

{ak-iXk-iXk + bk-iYk^iYk) + {akXkXk+i + bkYkYk+i) 



(01X1X2 + 611^1^2) + {o-2k-lX2k-lX2k + ^2fc-1^2fc-li^2fc) 
ll2k-lX2k-lX2k + b2k-lY2k-lY2k 
a.2kX2kX2k+l + b2kY2kY2k+l 



0-n-iXn-iXn + ^n-iYn-iYn 

Then starting from [Zj^^i — Zi^^Hq\, we can sequen- 
tially generate Z^-i — Zn, Zn-2 — Zn-i, Zk — 
Zk+i- Together with Zk, we can generate Z^+i, and 
hence Zk+i, ■ ■ ■ , Z2k-i, ■ ■ ■ , Zn- Then together with di^i + 
'^2*:-i-^2A;-i, we Can decouple and generate Zi. Analogously, 
starting from [Xk,Ho], we can generate Xi. Hence by Theo- 
rem]?] we have: 

Theorem 5. For an XXZ or XYZ chain of length N, with two 
local controls on the k-th spin. Hi = Z}^ and H2 = X^, and 
N ^ 2k + \ the system is controllable on the whole space. 

C. Controlling Z\ and X]^ 

Next we will see that two independent controls Zi and 
Xk on the first and the fcth spins are also enough for full 
controllability for a XYZ or XXZ chain, with Hamiltonians, 

N 

Ho = Y^ ajXjXj+i + bjYjYj+i + c^ZjZ^+i 

Hi = Zi, H2 = Xk. 

Earlier we showed that iJg and Hi can generate 

[Zi,Hq] -^Zj 

-^ajXjXj+i + bjYjYj+i 
-^ajXjYj+i - bjYjXj+i 

N 

-^^'^iZ.^Z.^^i = H3 
j 



and subsequently 

[Xk, H3] Ck-iZk-iYk + CkYkZk+i 
[Xk,a,k-iXk-iXk + bk^iYk^iYk] Yk-iZ^ 
[Xk,Yk-iZk] — > Yfe-iYfe 
[Yk-iYk, Ck-iZk-iYk + CkYkZk+i] Xk-i 

Continuing such process, we can sequentially generate Xk^2, 
Xk-z, ■ ■ ■ , X2, Xi. By Theorem]?] we have 

Theorem 6. For an XXZ or XYZ chain of length N, with two 
local controls Zi and X^, 1 < k < N, on the first and the 
k-th spins, the system is controllable on the whole space. 

V. Internal Symmetries and XX Chains 

In previous sections we have discussed the controllability 
results for XXZ or XYZ chains with one or two local controls. 
The fact that there are spin-spin interactions in all three 
directions is essential to guarantee subspace controllability or 
full controllability, and we expect that this will no longer be 
true if the spin chain is reduced to XX, XY or Ising types, as 
absence of the ZZ terms in Hq allows for more symmetries. 
Since the XX chain is particularly important as a theoretical 
model, we first investigate how the controllability properties 
change in this case. 

We start again with an XX chain and a single end control 
in the Z-direction 

N 

Ho^J2 i^r^Xn+l + r„y„+i ) ,Hi=Zi (7) 
n 

Theorem 7. For an XK chain of length N with a single 
local control on the end spin in Z-direction, the associated 
dynamical Lie algebra C is a faithful representation of u{N) 
in each excitation subspace Hk, k = 1, . . . , N — I. 

Proof In order to prove C is indeed isomorphic to 
u(A^), we determine all independent operators generated from 
the Hamiltonians, as we did in previous sections, and show 
that these operators satisfy the same commutation relations 
as the standard basis of u(A^) l^. However, to make the 
following analysis of faithful representations simpler and more 
convenient, we transform the original representation using the 
Jordan-Wigner(JW) transformation, a powerful tool initially 
developed in theoretical physics: defining 

am l[Z^{X- lY)^ 

where am are fermionic annihilation operators, with the canon- 
ical anticommutation relations 

the Hamiltonians are transformed into: 

N 

Hq^^ "fniaian+i + ajj+ia„) 

n=l 

Hi = a\ai 
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By calculating the commutation relations between the Hamil- 
tonians, we can verify the following identities: 

[iHi,iHQ] -> 2/12 := a|a2 - ajoi 
['iHi,yi2] xi2 := i{a\a2 + 4^1) 
[xi2,yi2] Z2 



ia\a2 



Then we can generate Hq = Hq — 71X12, which represents 
the system Hamiltonian for a chain of length — 1, and Z2 
amounts to the end control on Ho. Thus by induction we can 
sequentially generate 



yn,n+l • — "'n+l"-n: 
Xn,n+1 ■= i(a,t«n+l + ^n+l^'n), 

Zn ■■= i{alan), 



n=l,...,N -I 
n = l,...,N -1 
n = l,...N. 



Hence, we have generated three kinds of anti-Hermitian oper- 
ators: Xn,n+i, yn,n+i and z„, Satisfying 





~~ ymn 


(8a) 


[ymn 1 


— •^mn 


(8b) 




— Vinn 


(8c) 




— •^mn 


(8d) 






(8e) 



which are the same as the commutation relations satisfied by 
the standard basis of u(iV) |6|. Hence, we have C — u{N). 

Finally, in order to show £ is a faithful representation in 
each excitation subspace Hk, k — 1, • • • , A^ — 1, it is sufficient 
to show that the images of the generators x„,„+i, yn.n+i and 
Zn in Hk are nonzero. Choose a vector |a) E Hk in the 
computational basis with n and 71 + 1 positions as |0) and 
and define another basis vector G Hk such that /3 
only differs from a at n and n + 1 positions, with values |1) 
and |0). Restricted on Hk, we have 

{f3\yn,n+i\a) = {P\aian+i\a) = 1 
{I3\xn,„+i\a) = {I3\a'lan+i\a) = i 
= {P\ala„W) = 1 

Hence, on each Hk, k ~ 1, • • • , iV — 1, £ is a faithful 
representation of u(A^). ■ 
The above theorem implies that for XX chains with a Zi end 
control the system is controllable only on the single excitation 
subspace (or its mirror image, the iV — 1 excitation subspace). 
In order to derive full controllability we require more controls. 
If we add another local control Xi, we have the following 
result (proof in |5 1): 

Theorem 8. For an XX chain of length N with two local 
controls on the end spin, Hi = Z\ and H2 = Xi, the system is 
not controllable on the full Hilbert space and C = 5o(2N + 1). 

To better understand this result we find all operators in 
£ and the identify their symmetries. According to previous 
discussions, the system Hamiltonian of the XX chain, Hq, and 
Zi generate XkXk+i + YkYk+i, XkYk+i - YfeX^+i and Zk- 
On the other hand, Hq and Xi generate Z1Y2 and then Y1F2 
and XiX2- Similarly we can decouple and generate A'iy2 and 



YiX2- Also Xi and Zi will generate Yi, and together with 
Ho we can generate ^1X2. Continuing this process, we have 



[Z1X2, X2X:i + Y2Y:}] Z1Z2Y:} 

[Z1Z2Y3, Z3,] Z1Z2X3, 

[Zi ■ ■ ■ Zn ^2X^-1, Xn^iX]~^ + Y/v^iFat] — > 
[Zi - ■ ■ Zn^iYn, Zn] Zi - ■ ■ Zn-iXn 
\Z\ ■ ■ ■ Zj^iXj, Zi ■ ■ ■ Zk-iXk] XjZjj^i ■ ■ ■ 
■ Zk-iYk — ^ YjZj+i ■ ■ ■ Zk-iXk 



Z\ - ■ ■ Zn-\Yn 
Zk-\Xk 



XjZj+i 



Hence we can generate 4(^) operators of the form 
PjZj+i ■ ■ ■ Zk-iPk, where P can be either X or Y, and 
there are A^ operators Zi ■ ■ ■ Zk-iPk, where P = A' or 
Y. Together with A^ number of Zk, we have dim(£) = 
4(^) +2N + N = N{2N + 1) = dimso(2Ar + 1). This 
suggests that the Lie algebra is a representation of so(2A^ + l), 
although since the Lie algebra is not the maximal Lie algebra 
on the subspace, calculating the dimension is not sufficient. 
However, setting S = Y1X2Y3X4 • • • , it is easy to verify that 
we have HjS + SHj = for j = 0,1,2 and thus S is 
an internal symmetry. Depending on the number of F's in 
S, S is either symmetric or antisymmetric and thus defines 
an orthogonal or sympletic symmetry. This means that £ is a 
subalgebra of so(2^) or sp(2^~^), respectively, and therefore 
not controllable. Strictly, this is still not sufficient to conclude 
that C ~ so(2A^ + l). However, the table of irreducible simple 
subalgebras of su(2^) in [E] suggests that so(2A^+l) is indeed 
the only possible subalgebra that has the correct dimension and 
symmetry. 

Next, we study the XX chain subject to two controls, Zi 
and X2- In this case, it is straightforward to verify that 



S — X1Y2X3Y4 



(9) 



defines an internal symmetry. Depending on the number of Y's 
in iS", denoted by Ny, the symmetry can be either orthogonal 
or symepectic: if Ny is even, S is an orthogonal symmetry, 
and if Ny is odd, S* is a symplectic symmetry. 

Theorem 9. For an XX chain of length N subject to two local 
controls Hi = Z\, H2 = X2, the dynamical Lie algebra 



C 



so(2^) 
sp(2^-i) 



A^ e {4fc,4fc + 1 : fc > 1}, 
A^e {4fc + 2,4/j + 3 : /c > 0}. 



(10) 



By the definition of internal symmetry, any operator M E C 
also satisfies the symmetry: M'^ S + SM = 0. Define the 
symmetry-preserving set Q to be the set of all general Pauli 
operators M = P1P2 ■ ■ ■ Pn, where P = X,Y, Z, satisfying 
M^S + SM ^ 0. We have £ C ^ and span(^) C su(2^) 
forms a Lie algebra. In the following, we determine rank(t/) 
and show that all operators in Q can be generated from the 
Hamiltonians, i.e., £ — spa.n{Q). Denote by Nf. (A^o) the 
number of even (odd) positions of the XX chain; by Ue (rio) 
the number of A"'s and Y's contained in M at the even (odd) 
positions, and nz as the number of Z's contained in M. We 
have N = Ne + No, <ne < Ne, and < < A^o. 

Here is the rule for any M e Cy: if ng is odd, then nz 
is even; if rig is even, then nz is odd. For example, X1Z2, 
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and both preserve the symmetry, while Z1X2 and 

X1Y2Z3 do not. Notice that the X's and Y's contained in 
A4 are interchangeable at the same location to preserve the 
symmetry. For example, Y1X2X3 and Y1Y2X3 are also in Q. 
This fact is useful when evaluating rank(tj): for a given M E 
Q, containing p number of X's and F's, we can write down 
2P number independent operators in Q based on M. 

Specifically, for N ^ 2k, S = X1Y2X3Y4 ■ ■ ■ Xn-iYn and 
Ne = No = k. If JT-e = 0, then riz can be 1,3, ... ,k, and tIq 

can be chosen as (''^°), (^°). (^°)' Denoting Nn^=m as 

the number of M € Q with rig = m, Nn^=Q is equal to: 



20 



21 



3 / 



I 



Analogously, for rig 



1, nz can be 0,2,4,... and Nn^=i is: 



I 



(?)2-- [a°) + (^°)+--- + (t)J-(t)2 

Repeating such process, we find Nn^^m = (m°)2^^^2^°. 
Finally, for — k, if k is even, then should be odd, and 
Nn^=k is equal to: 



1 
1 

-,k-lf N, 



[( 



' \N-k- 



I)] 



"-JO 



)2 



a- 

(^■=)2^-i(2^ 



2^° - 1) 

Therefore, 

= 2^-12^° + + • • • + I 

= 22^-1 - 2^-1 - dim (so(2^)) 
If k is odd, then nz should be even, and Nn^=k is: 



{ k 



N-1 



2nt) 2"(?)[(V)+••• + (^."^^)] 



(? 



[( 



N-k- 




(N-k- 
\N-k- 



I)] 



+ --- + 2'=-Hf_",)Q+2nl°)0 



t)2- 



[ 



(? 



Therefore, 



rank(5) = ^ 7V„ 

m=0 



dim(sp(2^-i)) 



For iV = 2fc - 1, we have = k and ^ k - 1. 
Analogously, for both even and odd k cases, we can respec- 
tively find the value of rank(tj). Specifically, for odd k, and 



rank(g) = 2^^-^ + 2^-^ = dim (sp(2^-i)); for even k, 
rank(g) = 2^^-^^ - 2^-'^ = dim(so(2^)). Hence we have 
the following result: 

Lemma 1. For an XX chain of length N, with two local 
controls Hi = Zi, H2 = X2, when N ^ Ak or Ak + l, k > I, 
rank(g) = dim {so{2^)); when N = Ak + 2 or 4fc + 3, fc > 0, 
rank(g) = dim (sp(2^-i)). 

Next, in order to prove Theorem [9] it is sufficient to show 
that C = span(^). Then, since £ is a subalgebra of so(2^) 
or sp(2^-i) and dim(/:) = rank(^) = dim (so(2^)) or 
dim (sp(2^-i)), C must be so(2^) or sp(2^-i). 

[Ho,Zi]^ Zk,k^l,...,,N 

— >■ XkXk+i + YkYk+i,XkYk+i — YkXk+i 

[X1X2 + Y1Y2, X2] — 7- — > Y1Y2 — > X1X2 
[X1Y2 — Y1X2, X2] — > X1Z2 — > X1Y2 —7- Y1X2 

Thus, with the help of X2, we can decouple and generate P1Q2 
and similarly Q2P3, where P ^ X,Y andQ = X, Y, Z. Then, 

[X2X3, X3X4 + ^3^4] ^ ^2^3^4, [X2X3, X2Z3Y4] — >■ 

Y3Y4 X3X4 • • • Yn-iYn Xn-iXn- And simi- 
larly we can decouple and generate all XkYk+i and YkXk+i. 
Hence we can generate all P^P^-^, where P — X,Y. We can 
see that it is very important to first generate a single decoupled 
X1X2, and then using it we can decouple all the other PkPk+i- 
This is possible only when we have X2 initially. Hence, this 
availability of X2 control is essential for Theorem |9] Notice 
that as we can generate all Zk, for any M E C that contains X 
at the fcth position, i.e., M = ■ ■ ■ X^ • • • , we can also generate 
[M, Zk] — ^ Ms = ■ ■ - Yk ■ ■ ■ . Hence, it is sufficient to show 
the result for M only contain X and Z operators. From X1Z2: 

[XiZ2,Y2X3] — > X1X2X3, [X1X2X3, X2Y3,] — > X1Z3 
[XiZ3,Y3X4] — > X1X3X4, [X1X3X4, X3Y4] — > X1Z4 

XiZn 

[XiX2X3,YiX2] — > 2^1^31 [^1^3 5 ^l-'^3-'^4] ^ X1X4 

Therefore, we can generate all XiZ^, k = 2,...,N, and 
X1X4. Next, 

[X1X3X4, Z1X4] — ^ Y1X4 — ^ X1X4 

[XiZ4,XiX4] -^Y4^ X4 

This means that from X2, we can generate X4. Then we 
can imagine that the chain starts from the position 3, and we 
control Z3 and X4, and analogously we can generate Xq. This 
observation is very important, as it tell us whenever we can 
generate an operator with respect to positions 1 and 2, we 
can automatically write down the analogous operators we can 
generate at positions 3 and 4, or 5 and 6, etc. For example, 
for an XX chain with length N — 20, since we have shown 
that we can generate X1X4, we immediately know that we 
can also generate X^Xq, X^Xg, etc. 

Now we are ready to prove Theorem [9] by induction. For 
iV = 2, it is trivial to show C = span(C7) = sp(2); for iV = 3, 
we first list all elements of Q in Table |l] 
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1-body 




2-body 




3- body 


Z1Z2Z3, Z\X2Z-j., X1Z2X3, X1X2X3 



TABLE I: All symmetry preserving operators in Q foi N = 3. 

In Table |l] besides the operators we have already generated 
as above, we can also get: 

[Z1X3, Z2Y3] — Z1Z2Z3 [ZiZ22'3, I2] Z1X2Z3 

Thus, we have generated all symmetry-preserving operators 
for = 3, and hence C = span(tj) — sp{2'^). Similarly, for 
Af = 4, we can generate all symmetry-preserving operators in 
Q. In particular, for 2-body operators, we also need to show 
we can generate X3Z4, which is true, as 

[ZiA3,YiZ4] ^ A1A3Z4 
[A1A3Z4, A1Z3] — > Y3Z4 A3Z4 

For all other 3-body and 4-body operators, we can sequentially 
generate them from existing 2-body and 3-body operators. For 
example, we have 

AiA2A3,y3A4] — > XiX2Z^X4 — > X1Z2Z3X4 

[A1Z2Z3A4, Y1Z4] — > Z2Z3X4 

We can verify that all operators in Q can be generated for 
A^ = 4, and £ = ^=50(2^). 

Next, assuming that for all A^ < 2k, we can generate all 
operators in Q, we aim to show that for A^ = 2A; + 1 and 
N ^ 2k + 2, we can also generate Q. First, when A^ = 2fc + 
1, for any M = P1P2 • • • P2fc^'2fc+i € G, if P2fe+i is the 
identity, then M is also a symmetry-preserving operator for 
the chain with length 2k, and by induction assumption, M can 
be generated. If P2i:+i is not the identity, but Pi and P2 are, 
then as discussed earlier, M can be considered as a symmetry- 
preserving operator for the XX chain starting from positions 
3 to A^, with length 2k — 1, and by induction assumption, it 
can also be generated. The remaining case is when AI is an 
operator where both P2fc+i and Pi or P2 are not identity, then 
we can always choose a symmetry-preserving operator Mg = 
P1P2Q3Q4 from positions 1 to 4 such that Q3Q4 is different 
from P3P4. For example, if M = A1Z2A3A4 • • • , then we 
can choose Ms = X1Z2X3. Defining L = [M,Ms] ^ 0, 
1/ is a also symmetry-preserving operator for the XX chain 
from positions 3 to A^, and M can be generated from L and 
Ms, which can both be generated by the induction assumption. 
Thus, we can generate all operators in Q for A = 2A: + 1. For 
A^ = 2fc + 2, we repeat the same induction argument again 
and derive the same result. Thus, we have completed the whole 
proof of Theorem [9] 

Theorem |8] and Theorem |9] strongly suggest that for an 
XX chain, two independent local controls are not enough 
for full controllability. Moreover, we can numerically check 
for a given chain length A^ that two local controls with any 
configurations do not induce full controllability, i.e., for two- 
control problem, we will never find a full controllability result 
that is true for arbitrary A^. Nevertheless, it is worthwhile to 



notice that when A^ = 4fc + 2 or 4A: + 3 and £ = sp(2^"i), 
the system is pure-state controllable |12|, i.e., we can generate 
arbitrary target pure state from the initial pure state under 
these controls. Pure-state controllability can be useful for many 
practical applications. 

Finally, we find that three local controls are enough to make 
the XX chain fully controllable: 

Theorem 10. An XY or XX chain of length N is fully 
controllable under three local controls Zi, Xi and X2 on 
the first and the second spins. 

Proof: From previous discussion, Zi, Xi and Hq for 
XY or XX chain will generate Zj, j = 1,...,N, P1P2, 
P = X,Y,Z, and ajXjXj+i + bjYjYj^i. Together with 
Ai and A2 we can further generate ZiXj, ZiYj, and 
Z1Z2. We have [Z1A2, 02X2X3 - &2i"2X3] ^ ^1^2X3, 
[A1A2, 02^2X3 + 625^^3] — > X1Z2Y3, [Y1Z2, X1Z2Y3] — > 

^1^3' [-^1 ^2-^*^3 5 ^1^3] ~^ ■2^2-2^3' [y2y3:-^2] ■^^2^3' 

[Z2Y3,Z2Z3] X3. Hence, from Hq, Zi, Xi and A2, we 
can first generate P1P2, and then A3 and Z2Z3. Then together 
with Z3, we can generate all P2P3, P — X,Y, Z. Continuing 
such process, we can sequentially generate PjPj^i. Then the 
following procedure is exactly the same as the XYZ chain 
case with end controls Zi and Ai, and will generate all 
fc-body Pauli operators and the system is controllable, with 

/: = su(2^). ■ 

VI. Conclusion 

In this paper we have studied two kinds of Lie algebra 
symmetries, external and internal, and applied them to analyze 
controllability problems for spin chains. Specifically, for XXZ 
and XYZ chains under single end local control, the Hilbert 
space decomposes into a set of invariant subspaces and we 
have shown that the system is controllable on each of these 
invariant subspaces. These models are arguably the simplest 
type of quantum computers one can think of: a physical Hamil- 
tonian with a single switch giving rise to an exponentially large 
dynamical Lie algebra. We have also shown that this result 
stills holds when there is control leakage effecting neighboring 
spins, implying the robustness of this controllability result 
against control leakage which is common for practical control 
systems. 

We have also addressed the question of minimal control 
resources for full controllability on the entire space. For 
various cases, we find that two independent local controls are 
sufficient to make the system fully controllable for XXZ and 
XYZ chains. Finally, we have studied the effect of eliminating 
the ZZ-coupling that reduces the XXZ chain to an XX chain. 
We find that in this case for a single local end control, the 
dynamical Lie algebra is the same on each invariant subspace, 
and hence the system is only controllable in the first excitation 
subspace. Moreover, through the investigation of the internal 
symmetries, we have discussed the two-control problem of an 
XX chain, finding that two local controls are still not sufficient 
for full controllability in the XX-coupling case. However, we 
can still get pure-state controllability given two controls if the 
two controls are Zi and A2 and the chain length is A = Ak+2, 
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or 4fc + 3. In addition, we have shown that the minimal control 
resource for full controllability is three local controls acting 
on the first two spins. This analysis provides a better under- 
standing of the relationship between Hamiltonian symmetries 
and the system controllability, and minimal resources required 
to achieve certain levels of controllability, which is important 
for both quantum control theory as well as practical control of 
spin chains. Furthermore, the techniques developed to compute 
the dynamical Lie algebras explicitly by decomposing the Lie 
algebra into n-body interaction terms, which can be iteratively 
generated from the original Hamiltonians, will be useful 
to investigate controllability for other spin networks where 
full controllability cannot be inferred from the propagation 
property. 
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following we will only give the X^X^+i + YfcYfe+i terms. 
Next, we have 

[XkXkj^i + YkYk+i, Xk+iYk+2 — Ykj^iXk+2] 

-^XkZk+iXk+2 + YkZkj^iYk+2 

—^XkZk^iYk+2 ^ YkZk+iXk+2 

[XiY2 - YiX2,P{Z)] ^ {X1X2 + YiY2)Z3 
[X1Y2 - Y1X2, {X1X2 + YiY2)Z^] ^ {Zi - Z2)Z^ 
[X2Y3 — Y2X^, {Zi — ^2)^3] — > Zi{X2X:i + Y2Y3) 
[X2Y3-Y2Xs,P{Z)] 

-^Zi{X2Xs + Y2Y3) + {X2XS + Y2Y-i)Zi 

->(X2X3+y2i"3)^4 

[^2^3 - ^2^3, (^2^3 + Y2Y-i)Zi] ^ {Z2 - Z^)Z^ 

Continuing this process, we have 

[Xgn - >3^4, {Z2 - Z^)Z^\ ^ ^2(^3X4 + Fan) 

{X^Y^-Y^X^,R2\ 
^^2(^3X4 + Y^Y^) + (X3X4 + Y^Y^)Z^ 
^(^3X4 + Y^Y^)Z^ ^ (Z3 - Z4)^5 
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Appendix 

A. XX Z and XYZ Chain with Z\ Control 

We calculate the operators in £ generated by and H\ 
in (|3]i: 

[Zi,ffo]^^l>2-1^1^2 

\Z\,X\Y2 — Y1X2] X1X2 + Y1Y2 
[X1X2 + Y1Y2, X1Y2 — Y1X2] Z2 — Z\ ^ Z2 



Continuing the process, we can generate all Z-^, X^Xfe+i + 
YfeYfe+i, and P{Z) — J^f ^j'^jZjZj+i. Noticing that when- 
ever we can generate XkXk+i+YkYk+i we can also generate 
its conjugate term XkYk+i — YkXk+i, for simplicity, in the 



Thus, we can generate 

Zk-i{XkXk+i + YkYk+i) 
{XkXk+i + YkYi;+i)Zk+2 
Zk-i(Zk — Zk+i) 
{Zk — Zk+i)Zk+2 

Then we have 

[Zk-i{XkXk+i + YkYk+i),Xk-iYk — Yk-iXk] 
-^{Xk-iXk+i+Yk-iYk+i) 

[Xk-iXk + Yk-iYk, XkYk+1 ~ YkXk+i] 
-^{Xk-iXk+i+Yk-iYk+i)Zk 

[{Xk-iXk+i + Yfe_iYfe+i)Zfe, {Xk-iYk+i — Yk-iXk+i)] 
— >(^fe-i — Zk+i)Zk = Zk-\Zk — ZkZk+i, 

and together with P{Z), we can decouple and generate all 
ZkZk+i- Moreover, we have 

[XkYk+i — YkXk+i, Zk+iZk+2] 
-^XkXk+iZk+2 + YkYk+iZk+2 

[XkXk+iZk+2 + YkYk+iZk+2-, Xk+iYk+2 ~ Yk+iXk+2] 
-^XkXk+2 + YkYk+2 

[XkYk+2 ~ YkXk+2, Zk+2Zk+3] 
-^XkXk+2Zk+3 + YkYk+2Zk+3 ^ XkXk+3 + YkYk+3 

Continuing this processing, we can generate all XjXk + YjYk- 
Furthmore, from XjXk + YjYk and X^Xk + Y^Yk we can 
generate {XjXm^YjYm)Zk- Hence we can generate all {Zj — 
Zm)Zk, and then any decoupled term ZmZk- 

Thus we can find all fc-body operators: first, Zk is the 
only 1-body operator we can generate. Then we can generate 
XjXk + YjYk and XjYk — YjXk as the 2-body operates. 
Next we can generate the 3-body XYZ-mixed operators 
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{XjXm + YjY„i)Zk, from which we the 2-body Z operator 
{Zj — Zm)Zk- Then based on the existing operators, we can 
generate the following 4-body XYZ-mixed operators and 3- 
body Z operators: 

[{XjXk + YjYk)Zm, XmYn — Y,nXn] 

^{X,Xk + Y,Yk){X„,Xn + Y^Y,,) 

[{XjXk + YjYk)Zm, ZkZn] — > {XjYk — YjXk)Z,nZn 

[XjXk + YjYk, {XjYk - Y,Xk)Z^Z^] ^ [Z, - Zt) 

Continuing such process, we can sequentially generate t- 
body XYZ-mixed operators and (£ — l)-body Z operators. 
When i is even, we can generate: 

(^mi^m2 ^~ YjY^-^YjY^^^Zrf^i^ • • • Zjji^ 

(^mi^m2 Yjji-^Yjyi^'j(^Xjyi^Xjyi^ -\- Y^i^Y^i^^ Z^i^ • • • Zjyi^ 



Since the 2N number of variables Xj and i/j are linearly 
independent, the rank of E over M is the same its rank over 
C, i.e., ranki{(£') = rankc(i?). Next, over the field C, we 
can apply the following reversible transformations: zj = Xj + 
iyj and z* = Xj — iyj, and then all elements in E can be 
expressed as polynomials over C in terms of Zj and zj!. This 
is equivalent to considering raising/lowering operators in the 
algebra. Specifically, 

XjXk + VjVk = Re(z*Zfe) 
^iVk - UjXk = lin{z*Zk) 



When i is odd, we can generate: 

(^mi^m2 ^~ Yf^i-^Yyy^^^Zj^-i^ • • • Z^^ 

(^mi^m2 Yjyi-^^Yjyi^^(^Xjyi.^X^^ -\- Yyyi^Y^fi^^Zyy^^ ' ' ' Z ^ 



Then any element q{mi,m2) ■ ■ ■ q{m2p-i,ni2p) in E can 
be rewritten as Qfe^^ma) ' ' ' <3«2p-i^™2p)' where Q{z) 
is an operation that takes either the real or the imaginary 
part of z. Hence, all elements in E can be rewritten in 
terms of 2N number of independent complex variables: 
zi, . . . , zjv and z'l, . . . , z'^. Next, we can show that the 
space generated by the set E is the same as the one 
generated by the set F whose elements are in the following 
form: z;;^^™^ • • • z;;,^^_^z„,,^, where ruk G [1, . . . , N], 
k = l,...,2p. In order to see this, we will show that 
^7711^^2) ' ' ' Qi^rn2p-i^^2p) ^ can be generated by the 
elements in F: from the sum of or the difference between 



generate Q{z*^^z^^)z*^^z^^- ■ ■ z. 



■ z" Zm^ , we can 



Continuing 

{X„,,X„,2+Y,n^Y^^)---{X,r,,-2Xmt~i+Y^^_2Ym,-i)Z,nt^^c\\ proccss, we can use polynomials in the form of 



'3(^mi^ni2)^m3^m4 ' ' ' 



to generate elements 



B. Rank of M 

Theorem 11. For given N and p with N > 2p, the rank of 
the p-pair operator set E^^p is equal to (^) (^~^)- 

Since, when evaluating rank(i?7v,p), only linear relations 
between the operators of Em,p are involved, we can consider 
every element of Em,p as a polynomial in terms of 2N 
variables, and transform the original problem into evaluating 
the rank of a set of polynomials. Specifically, for positive 
integers N and p with N > 2p, let E be the set of any 
polynomials in terms of 2N variables xi,X2, ■ ■ ■ ,xn and 
2/1, 1/2, • ■ ■ , UN, satisfying the following form: 

P[xi,X2, ■ ■ .,XN;yi,y2, ■■■^Vn] 
=g(mi,TO2)(7(m3,TO4) • • ■ q{m2p-i,m2p)q{m2p-i,m2p) 

where k) can take two forms, either q{j, k) — XjXk + yjyk 
or q{j, k) = xjyk — yjyk, and TOfc's are distinct from each 
other, with G [1, ■ • ■ , N], k = 1, . . . ,2p. In other words, 
any element in i? is a product of p terms, each taking the form 
XX + yy or xy — yx. For example, the following polynomials 
are in E: 



in the form of Q(4ii^™2)Q(^m3^m4) ' ' ' 4i2p-i^m2p- 
Continuing such process we can finally generate 



On the other hand, reversing 



^m2p-l^'"2p 



Therefore, we 



such process, we can generate z^ 

from g(z;;,^ 

^m2p-i ^^2p ) 

have shown that span(i?) = span(F) over C, inducing 
rankc(£^) = rankc(F). 

Next, we evaluate the rank of F over C. Since the 2N 
variables zj and z* are linearly independent, all elements in 
F, as in the product form of z*^^ Zm^ • • • z,^^^^ z,„,p, are hence 
independent as well. In order to obtain an element in F, we 
choose p number of z^'s from the indices j e {!,..., N} and 
choose p number of z* 's from the remaining N — p number 
of indices. Then the number of elements in F is 
which equals to rankc(-F) ~ rankc(i?) = ranki{(-fc). Thus, 
we have proved rankR(i?)(^) (^jT^)' '^^'^ hence the result of 
Theorem [TT| follows. 



{xiX2 + yiy2){xzXi + 2/32/4) ■ 
(a;i2;2 - yiX2){x'i,yi - yax^) ■ 



ix2p-lX2p + y2p-iy2p) 
{X2p^iy2p - y2p-lX2p) 



C. Dimension of C 

After we derive all the fc-body operators in £, we evaluate 
dim(£), where the following identity is involved: 



Lemma 2. 



As discussed earlier, the total number of polynomials in 
E is p!(^)(^p^^)- However, not all of them are linearly 
independent, and we aim to evaluate rank(£') over JR. 



2N\ L^^J N\ 



E 



Nl ^^pP{N-2py. 



2N-2p 



Proof: Define a polynomial f{x) = {x + 1) . Then the 



term has coeffient { We have 

f{x) = (a;2 + 2x + if 

^-^ v\q\r\ 

^ jjIg'ICiV-p-g)! 



Hence, 



p!(Ar- 2p)!p! 

L^/'J /VI 
^=0 l'!'(^-2p)! 



